Which Quantum Evolutions Can Be Reversed by Local Unitary Operations? 
Algebraic Classification and Gradient-Flow-Based Numerical Checks 
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Generalising in the sense of Harm's spin echo, we completely characterise those unitary propagators 
of effective multi-qubit interactions that can be inverted solely by local unitary operations on n qubits 
(spins-^). The subset of U £ SU(2 n ) satisfying U~ = K\ UK2 with pairs of local unitaries K\ , K2 £ 
SU(2)® n comprises two classes: in type-I, K\ and K2 are inverse to one another, while in type-II 
they are not. Type-I consists of one-parameter groups that can jointly be inverted for all times t 6 R 
because their Hamiltonian generators satisfy KHK^ 1 — Adic(J?) = —H. As all the Hamiltonians 
generating locally invertible unitaries of type-I are spanned by the eigenspace associated to the 
eigenvalue —1 of the local conjugation map Adx, this eigenspace can be given in closed algebraic 
form. The relation to the root space decomposition of sl(N, C ) is pointed out. Special cases of type-I 
invertible Hamiltonians are of p-quantum order and are analysed by the transformation properties of 
spherical tensors of order p. Effective multi-qubit interaction Hamiltonians are characterised via the 
graphs of their coupling topology. Type-II consists of pointwise locally invertible propagators, part 
of which can be classified according to the symmetries of their matrix representations. Moreover, 
we show gradient flows for numerically solving the decision problem whether a propagator is type-I 
or type-II invertible or not by driving the least-squares distance Ki — e +lti/ ||2 to zero. 



PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Yz, 03.67.Pp; 33.25.+k, 76.60.-k; 82.56.-b 
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Introduction 

Richard Feynman's seminal conjecture that quantum 
systems may be used to efficiently compute and predict 
the behaviour of other quantum systems Q has inaugu- 
rated branches of research inter alia dedicated to Hamil- 
tonian Simulation [E IE EL IE IE • Actualhy, backed by 
the considerations by Manin J8j, Bennett [3] and others, 
it initiated efforts to explore the power of quantum com- 
puting. Soon thereafter, sets of universal one- and two- 
qubit gates were found [l(j which allow for decomposing 
any unitary representation of a quantum computational 
gate into elementary universal ones. 

Exploring computational complexity as well as devis- 
ing timcoptimal realisations of given quantum algorithms 
by admissible controls has therefore become an issue of 
considerable practical interest, see e.g. ^301- I n par- 
ticular the number of computational steps required to 
implement a quantum gate or to simulate its Hamilto- 
nian is a measure of the actual cost to put the gate or the 
simulation into practice. A specific question is, whether 
the sign-inverted Hamiltonian —H can be simulated with 
only +H and a set of given control Hamiltonians being 
at hand. The work of Beth et al. has addressed this 
problem for pair interactions to give bounds on the time- 
overhead HE Hi EE] required for doing so. In view of ef- 
fective multi-qubit interactions, here we go beyond pair 
interactions and classify those Hamiltonians that allow 
for simulating — H by H and local controls with exact 



time-overhead 1. 

This is of practical relevance, because when simulating 
quantum systems one often faces two-part generic tasks: 
(i) let certain interactions evolve while (ii) other effec- 
tive multi-qubit interactions shall be suppressed. The 
latter may be achieved by decoupling, but often it suf- 
fices that unwanted interactions cancel at a certain time, 
e.g. right at the end of an experiment, which is to say 
they are to be refocussed by inverting them at suitable 
points in time. This is important for instance to avoid 
undesired or dissipative coupling of a quantum system to 
its environment or bath |l6j | . Man y te chniques have been 
developed in magnetic resonance [jj], HE] on the basis of 
averagc-Hamiltonian theory [l9l ]. Moreover, local inver- 
sions arise in the context of LOCC, i.e. local operations 
and classical communication pfl ]. 

Let the operator T(t) denote time translation by t 
while Q represents time reversal. Following Wigner [211 
|23 | in these very general terms, one immediately finds 

T(t)o6oT(t)o6 = l«9oT(t)oe = r 1 (t) =T(-t). 

(1) 

Now imagine time translation is accompanied by the 
Hamiltonian unitary evolution of some quantum interac- 
tion H according to U(t) = e~ itH for all t. Then Eqn.ffl 
turns into 



e o u(t) o e = u(-t) 



(2) 
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Clearly time reversal itself is an unphysical operation, 
however, there are manipulations that bring about ef- 
fective time reversal for evolutions of certain quantum 
interactions, the most prominent early example of which 
being Hahn's spin echo |23| . 
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K ■ U(t) ■ K 1 = U(-t) = U(ty l K x ■ U(x) ■ K 2 = U(-x) = U(xy l 

Figure 1: There are two non-trivial instances of locally invertible quantum evolutions: (a) in Type-I the interaction can be 
refocussed by local unitary conjugation with K £ SU(2") for all times tel. (b) Type-II locally invertible interactions can 
only be refocussed at specific times r £ E since K\ and K2 are not inverse to one another. 



Generalising the sense of Hahn's spin echo, here we 
ask for which (non-trivial) Hamiltonian evolutions effec- 
tive time reversal can be obtained by local unitary op- 
erations, in other words which Hamiltonian evolutions 
can be entirely refocussed by framing them solely with 
local unitaries. As illustrated in Fig. ^ for n qubits, by 
this we mean: a unitary quantum propagator or gate 
U := er ltH £ SU(2") (with non-zero t £ M) is invertible 
exclusively by local unitary operations, if 

3K U K 2 £ SU(2)®" : K ie - %tH K 2 = e +itH . (3) 

Whether or not K\ and K 2 are inverse to one another 
has implications for universality, as will be shown next. 



Case Distinction 

As illustrated in Fig.^ locally invertible propagators ex- 
ist in two types: 

Lemma Either e~ ttH is trivial and self-inverse, or 
(i) it is type-I invertible in the sense 3K £ SU(2)®" : 
KHK" 1 = -H so Ke- itH K- x = e +itH jointly for all 
t £ M , or (ii) it is type-II invertible such that at some 

HtH with 



(but not all) points r in time K\e lrH K 2 
K U K 2 £ SU(2)® ra and K 2 £ K~ 

Since self-inverse cases with U 2 = 1 (as in the quantum 
computational CNOT, swap and Toffoli gates) are triv- 
ial from the point of view of inversion they will be not 
discussed any further. In order not to bother the reader 
with a stumbling block, the case distinction of Lemma 
will be proven in Appendix A. 



Organisation and Notation 

Rather, for the sake of illustration, we first address 
the problem of type-I invertibility from a geometric point 



of view leading to Lie-algebraic terms. A superoperator 
representation of the adjoint mapping then turns out to 
reduce the problem to a simple cigcnoperator calcula- 
tion that can be solved algebraically in closed form. This 
paves the way to discuss pair interaction Hamiltonians. 
Connecting several pairs then leads to assess local invert- 
ibilty in terms of coupling graphs in the case of multi- 
qubit systems. The transformation properties of multi- 
qubit interactions seen as spherical tensors of different 
quantum orders p allow for treating the problem in its 
normal form, namely invertibility by joint or individual 
local ^-rotations. Moreover, the quantum order p relates 
to the roots of the standard root-space decomposition, 
which gives necessary and sufficient conditions for local 
invertibility. After relating the findings to Cartan-like 
decompositons induced by the concurrence Cartan invo- 
lution linked to time-reversal symmetry, we finally es- 
tablish the relation to global minima of the least-squares 
distance over the restricted group of local unitaries. So 
invoking the norm property associated with the Frobe- 
nius distance, e~ ltH is locally invertible if and only if 



mm 

K 1 ,K 2 eSU(2) 



\Kie 



-itH 



Kn - e 



+itH\\2 _ 







(4) 



which can readily be decided on numerically by a gradient 
flow restricited to the local unitary group. 

Type-II invertibility is first treated by establishing 



Table I: Types of locally invertible interaction evolutions 



local invertibility 


K 2 = K- 1 


K 2 + Kl 1 


jointly for all t £ R 


Type-I 


{} 


point wise for some r £ K 


self-inverse 


Type-II 
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symmetry properties for the matrix representations of 
the interaction Hamiltonians. Then a system of two cou- 
pled gradient flows is devised to solve the problem nu- 
merically. It can be seen as a flow for the singular-value 
decomposition (SVD), yet restricted to local unitaries U 
and V. 

Throughout the paper, we use the following notation: 
G := SU(2"), K := SU(2) <g> SU(2) ® • • ■ ® SU(2) =: 
SU(2)®" for the Lie groups of the special unitaries and 
local unitaries as well as g and t for their respective Lie 
algebras. Elements of G, K, g and 6 arc written as 
G,K;g,k with the only obvious exception of expressing 
Hamiltonians by iH G g. 



I. LOCALLY INVERTIBLE ONE-PARAMETER 
UNITARY GROUPS 

Prelude: Geometry 

In the single-qubit case, a spin- 5 rotation U(n,<j>) by 
some angle <f> about the axis n is — most intuitively — 
inverted by a 7r-rotation about some axis orthogonal 
to n according to 



Uin,^)- 1 = U(n ± ,n) U(n,4>) {/(n^.Tr)" 



(5) 



Let the Lie algebra su(2 n ) be spanned by some orthonor- 
mal basis set {a,j}. Then the generalisation of rotations 
to higher dimensions, e.g. n qubits (spins-^) is straight- 
forward: replace the rotation axis by the subspace of the 
Lie algebra su(2" ) that is invariant under the action of 
H 

I H := span R {a, e su(2") | [ aj ,H]=0] , (6) 
and consider its orthocomplement in su(2") 

ijj := span R { aj G su(2™)| a-j & I H } ■ (7) 



Making use of the Hilbcrt space structure [24| , every H 
induces a specific decomposition 



su(2 n ) 



I'M 



(8) 



This setting already implies a particularly simple and 
illustrative first characterisation of locally invertible uni- 
taries, which, however, is not yet complete: 



Lemma 1 For a propagator U 



itH 



to be locally 



invertible for all t, the orthocomplement 1^ to its in- 
variant subspace in su(2 n ) necessarily has to comprise at 
least one local effective Hamiltonian k with K := e 
SU(2) 8n . 



— ik 



Proof: Assume there were no local Hamiltonian k G ig : 
then all the local unitaries K = SU(2)®" would be an in- 
variant subgroup to SU(2") under the action of the one- 
parameter unitary group {U = e~ ltH \ t G R }. In turn, 
there would be no local unitary to invert a propagator 
U = e- itH for all t. ■ 



Lemma 2 For a U = e 



-itH 



to be locally invertible at all 



t, it is sufficient that there is a local Hamiltonian k in 
the orthocomplement Jjj so that the double commutator 
of H with k reproduces H, i.e. [k, [k,H]] = H. 

Proof: In the first place, note that 

span K {£f, k, i[H, k] } = su(2) implies [k, [k,H]] = H, 
whereas the converse does not necessarily hold. How- 
ever, the condition [k, [k, H]] = H suffices to define an 
analytic function 



with the derivatives 



d 2 f 



icj>k. 



[k,H] 



-icpk 



[k,[k,H]]e 



i(pk 



(9) 
(10) 

(11) 



H and M- 1 , „ 

d<j> I 0=0 



The boundary conditions /(0) 
—i[k, H] allow for expressing the function f{4>) as 



H coscj) — i[k, H] sine 



For. 



7r and K 



-ink 



G SU(2)®" one finds H 1 



(12) 
-H 



so e is locally inverted. 



However, for obtaining a both necessary and sufficient 
condition, it seems one has to sacrifice the illustrative 
simplicity of geometry. 

Lemma 3 For a U = e~ ltH to be locally invertible at all 
t, it is both necessary and sufficient that there is a local 
Hamiltonian k in the orthocomplement 1^ and a suitable 
4> G [0, 2?r[ with 



OO ^ 

£7 



1 1 ad 



-H 



(13) 



Proof: Immediate consequence of the well-known iden- 
tity 



e x Ye- x = 



OO -j 



(14) 



e=o 



where ad^-(F) is the £-fold commutator of X with Y, i.e. 
[X, [X,... [X, Y] ■ ■ ■ ]] and ^ad£ := 1. Note this includes 
Lemma 2 as a special case. I 



Algebra I: Eigenoperators 

To begin with, recall the well-known fact that for two 
matrices A, B to be similar i.e. XAX" 1 = B with a 
non-singular X, their eigenvalues have to coincide. Thus 
a Hamiltonian H is invertible by unitary conjugation, if 
and only if its non-zero eigenvalues (including multiplic- 
ity) all occur in pairs of positive and negative sign. This 
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is a necessary and sufficient condition for inversion under 
some U G SU(2"), whereas for local unitary inversion 
by a K G SU(2)®" being a special case it is merely a 
necessary one. 

Complete Basis for Locally Invertible Hamiltonians 

Although the decomposition of the algebra su(2)®" into 
invariant subspace and orthocomplement is illustrative, 
it is very tedious to be carried out case-by-case for each 
and every given Hamiltonian H. Rather, in order to ob- 
tain constructive parameters, we will turn to the group 
SU(2)®" of local unitaries and give a basis set to its 
cigenspace in which all locally invertible Hamiltonians 
(of type-I) can be spanned. To this end, observe that 
due to the series expansion, 



Ke- itH K- x = e~^ KHK " 
KHR- 1 = Ad K (H) = -H 



e +itH V t g M (15) 



(16) 



where the latter in turn is equivalent to the series expan- 
sion of Adx in Lemma [3] and Eqn^] Moreover, by the 
Kronecker product and the notation of a matrix as a vec- 
tor ('vec') consisting of the matrix columns stacked one 
upon another l2q | one has illustrated the following 
obvious necessary and sufficient criterion for local invert- 
ibility given as assertion (l)in the following 

Lemma 4 (1) The propagator e~ ltH is locally invertible 
for all t G K if and only if vec H is eigenvector of AdK 
(here represented as K* ®K) to the eigenvalue — 1: 



(K* (g> K) vecH = -vecH 



(17) 



(2) The eigenspace to the eigenvalue —1 spans all the 
locally invertible Hamiltonians, and it can be given in 
closed algebraic form by recursively making use of the 
eigenvectors in SU(2), as K G SU(2)®". 

Proof: for assertion (2), we give a constructive proof in 
view of explicit applications. 

Eigenvectors of the AdK-Superoperator K* (g) K: 

For any local unitary K G SU(2)® n , the superoperator 
K* ® K is just a 2n-fold tensor product of unitary 2x2 
matrices. Using the quaternion parameterisation 



U := cos f 1-i sin ^(n x a x +n y a y + n z a z ) G SU(2) (18) 

+ •§. 

with J2v=x y z n l = 1> tnc eigenvalues A± = e 1 2 (let 
j3 ^ 0) are associated with the orthonormal eigenvectors 

1 



v+ 



where the limit n. 
and u_ = ( J ) . 



^2(1 + n z 
1 



1 

1 4 



n z 
- in. 



(19) 
(20) 



— ^ — 1 is uncritical: one finds v + = ( J ) 



One Spin-^ Qubit 

The AdK-superopcrator (A"* ® K) for a single spin 
qubit thus shows the four eigenvalues A±A ± = e TZ 2 e ±l 2 
(being either 1 or e* 1 ' 3 ) associated with the four orthog- 
onal eigenvectors v± <S> v± . Consequently the eigenspace 
to the overall eigenvalue e ±l13 — — 1 is spanned by the 
basis set 

E { ~ } := {v*_ ® ® (21) 

while the eigenbasis to the overall eigenvalue +1 reads 

:= {u* <X> ® «_} . (22) 

Remark 1 For fixed parameters (n Xl n y ,n z ) one finds 
however, that every element in su(2) 
can 6e spanned in both E^ and E^ + ' : e.g., vec(a z ) may 
be expanded in E^' by (n X: n y ,n z ) — (cos#, sin#, 0) ; 
while in E^ + > the expansion requires (n x ,n y ,n z ) = 
(0,0,1). This re-expresses the trivial fact that a z is in- 
verted by any n rotation about some axis in the xy -plane, 
whereas it is invariant under z rotation. 

For completeness, in the limit (3 — > define E^ := 
^= {vec a x , vec a y , vec a z } . 

Two Spin-^ Qubits 

For two qubits, the eigenbasis to the overall eigenvalue 
\\±\2± A 1± A 2± = —1 consists of vectors of the following 
subtypes: 

Subtype Embedding of the two limiting 1-spin cases 
with \Pi \ = 7r or |/32 1 = tt 

:= {e[~ ] ® vec(l) 2 } U {vec(l) 1 ® E [ f ] } 

Subtype 1 Inversion of one spin or the other spin with 
|/3i|=7r,/3 2 = 0or /3i = 0, |/3 2 | = tt 

E„, i={Et ) ®4° ) }u{4 0) ®4~ ) } 

Subtype 2 Rotation on both spins with \/3i \ — I— | /?2 1 = 
7r(mod 27r) and fli,02 

Eflufh :={4 _) ®^ _) } 

Subtype 3 Rotation on one spin, commutation with 
the other spin where |/3i| = tt, f$i ^ arbitrary, or 
I $2 1 = 7r, /3i 7^ arbitrary 



n Spin-^ Qubits 

The generalisation to n qubits with 



■1 (23) 



is obvious, because the construction follows the pattern 
described by the indices to the cigenspaces. One may go 
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Table II: Type-I Invertibility of Elementary Pair Interactions 



Pair 
Interaction 


Expansion 
in Subtype 


Type-I Local 
Inversion by 


Symmetry 
Class 


ZZ 




tt(_L 1) or: tt(1 _L) 




XX 


E Pl,02 
E TT,0l\ 


n(zl) or: ix{lz) 

-/&(!*) 
[e.g.: f(2l-U)] 
tt(_L 1) - tt(1 H) 




XY 


-E-ir.O 


7r(zl) or: n(lz) 
n(xl) ± 7r(ly) 
or: n(yl) ± 7r(lx) 




X(-X) 


E n,0 
E Pl,02 

E TT,!3\\ 


7r(zl) or: n(lz) 
/3 1 ( 2 1)+/3 2 (1 2 ) 
[esp.: §(zl + L?)] 
tt(_L 1) +tt(1 h) 


• • 


XXX 
XXY 
XYZ 


none 
none 
none 







in Tab. [H] while the Ising-ZZ interaction can only be 
expanded in the eigenspaces of Subtype 1, i.e. £V,0j 
Heisenberg-XX and XY interactions allow for expan- 
sions in Subtype 1 as well as Subtype 2 (Ep u p 2 ~). 

For brevity, in the table we use the short-hand nota- 
tion (zl) for i(£7 z <g>l), and (zl±lz) for \{a z ®t±t®a z ). 
as well as (_L 1) for \{a x cos<p + a y sin^>) <g> 1 and anal- 
ogously with reference to some fixed </> we write (H) and 
(h) for \{o x cos(4> ± ?) + Uy sin((/) ± £)). For example, 
the Heiscnbcrg XX interaction can of course be inverted 
by 7r z-pulses on one or the other qubit, but also by an 
antisymmetric z-rotation on both qubits, where the ro- 
tation angle is ft on qubit 1 and ft — it on qubit 2. Note 
that inverting generic ZZ, XX, and XY interactions 
requires pulses that are non- symmetric with regard to 
permuting qubits 1 and 2. In view of convenient exten- 
sions to networks of pair interactions, we write • — o for 
a pair interaction of two qubits that is inverted by such 
non-symmetric local pulses. The only exception of dif- 
ferent symmetry is the Heiscnbcrg X(— X) interaction, 
since it can also be inverted by a permutation symmetric 
^ z-pulse on both of the qubits expressed by • — • . 

Note that none of the Heisenberg XXX or XXY or 
XYZ interactions is typc-I invertiblc by local unitaries, 
because their interaction Hamiltonians already fail the 
simple necessary condition of being invertible over the 
entire unitary group: their non-zero eigenvalues do not 
occur in pairs of opposite sign. For instance, the eigen- 
values to the XYZ interaction Hamiltonian Hxyz '■= 



a{<j x 



; ) + P(o- y <g> o-y) + 7(er 2 ® a z ) read 



from n — 1 spins to n spins by adding the nth index from 
the set {1, 0, /3, /3\\} to each of the previous n — 1-spin 
cases according to the subtype of embedding. Subtype 
means expand E n _\ to £7„_i ® veel; subtype 1 gives 
£ n _i <g> £K°); subtype 2 leads to £ n _i <g> E^; subtype 3 
results in E n _i ® E^+K 

In view of constructive results, the above subtypes have 
been expressed in terms of sets of consistent rotation pa- 
rameters (n x ,n y ,n z )i and rotation angles on every 
spin t. A locally invertible Hamiltonian has to be ex- 
pandiblc in at least one set of these self-consistent pa- 
rameter sets. ■ 

In larger spin qubit systems, these checks may become 
increasingly tedious. However, physical problems are of- 
ten confined to special settings: a Hamiltonian may be 
constituted by pair interactions, or in other instances, 
a Hamiltonian may be made up of terms that can be 
grouped in combinations of interactions transforming like 
spherical tensors of various p-quantum order. For these 
two practically relevant cases, we present more conve- 
nient methods. 



Ising and Heisenberg Pair Interactions 



Ai = 
A 2 = 
A 3 = 
A 4 = 



+a + j3 - 7 
-a + /3 + 7 
+a - p + 7 
—a — ft — 7 



(24) 



with a, ft, 7 € R . Clearly, unless at least one of the pa- 
rameters {a, ft, 7} vanishes, there are no pairs of opposite 
sign thus limiting the type-I invertible interactions to ZZ 
or XX or XY type. 



Coupling Graphs for Networks of Pair Interactions 

Coupling networks made up by pair interactions be- 
tween qubits can conveniently be represented by graphs: 
each vertex denotes a qubit, and an edge connecting two 
qubits k and / then corresponds to a non vanishing pair 
interaction or coupling Jki- These may take the form 
of any of Ising or Heisenberg type interactions described 
before. We will discuss connected graphs that do not 
necessarily have to be complete. 

As will be seen next, interactions with coupling topolo- 
gies of bipartite graphs have special properties. 



The pair interactions of Ising and Heisenberg type can 
easily be related to the —1 eigenspaces as summerised 



Lemma 5 (Variant to Beth (T^) The evolution un- 
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Figure 2: In a bipartite coupling graph, only vertices of dif- 
ferent colour o or • are pairwise connected. 



and analogously write Fz or Fz if the sum just ex- 
tends over all qubits coloured o or •, respectively. 

(i) For k £ ] — 1; +1] a bipartite coupling topology suf- 
fices to allow for the inversion Hxy i— ► —Hxy by the 
rotations nF^ or ttFz*\ A bipartite topology is also 
necessary, since in general no permutation-symmetric in- 
version of -Hxy exists (see Tab. ITT1). 

(ii) Cleary, also in the special case k = — 1 a bipartite 
coupling graph suffices. However, it is not necessary, be- 
cause a ^ z-rotation on all the qubits (^ F z ) is invariant 
under qubit permutation and thus does the same job on 
any connected coupling graph without requiring the dis- 
tinction of a bipartite topology (cp the permutation sym- 
metric inversion of the X(—X) interaction in Tab. lTT)l . ■ 



der Ising Z Z -interactions 

n 

Hzz ■= J kl \okz ® <Jlz (25) 

k<l 

is type-I invertible by local unitary operations if and only 
if its coupling topology of non-vanishing couplings Jki 
forms a bipartite graph. 

Proof: 

(i) For Hzz l_ * —Hzz it is sufficient that each of the 
edges of the coupling graph is inverted. Using local ac- 
tions on the vertices, this means every vertex of either 
the type • (or o) has to be inverted an odd number of 
times, while the other type o (or •) remains invariant i.e. 
is inverted an even number of times incl. zero. 

(ii) Not only is this condition sufficient, it is also nec- 
essary: assume there were edges connecting two vertices 
of the same type (either • or o). Then the couplings 
depicted by such edges would not be inverted, as they 
flip their signs twice (or an even number of times) thus 
remaining effectively invariant. ■ 



Lemma 6 The evolution under the Heisenberg 
XY -interaction 

Hxy ■= 7T ^ J kl \ (o-fex ®o-ix + k o- ky ® o ly ) , (26) 

k<l 

where n € [— 1;+1] is type-I invertible by local unitary 
operations if and only if (i) either its topology of non- 
vanishing couplings Jki forms a bipartite graph or (ii) 
k = —1, in which case the coupling topology may take the 
form of any connected graph. 

Proof: 

Let F v with v £ {x,y,z} denote the sum over n qubits 
with the Pauli matrix o\P in the £ th place 

n 

(27) 



Examples of Pair Interaction Hamiltonians 

For instance, neither Ising ZZ-coupling nor the Heisen- 
berg XX and XY interactions on a cyclic three-qubit cou- 
pling topology (C3) are type-I invertible, because C3 is 
clearly not bipartite. However, also on C3, the Heisen- 
berg X(-X)-interaction is type-I invertible, as will be il- 
lustrated below in the section on gradient flows. 



Extension to Effective Multi-Qubit Interactions 

In multi-qubit effective interaction Hamiltonians on a 
coupling graph G, the interaction order I (e.g. I = 3 for 
H c ff = cr z <E> o- z <E) o\z/2) may be used to group terms of 
different order. To each order, there is a subgraph G me . 



Lemma 7 Let H c ff be an effective multi-qubit interaction 
Hamiltonian constituted by the t-interaction terms on the 
me. subgraphs. 

H e g = Hi. me (28) 

where £ runs over the interaction orders and mi com- 
prises all the l-order interaction terms on subgraphs G mi ■ 
Then H c ff is locally invertible of type-I if and only if 
its constituents on the G me are all simultaneous eigenop- 
erators of some Ada to the eigenvalue — 1. 

Proof: 

The interaction order as well as the assignment to the 
subgraphs is Adx invariant. ■ 
In simpler cases, the grouping may help to find local 
inversions by paper and pen. However, more complicated 
multi-qubit interactions can be treated by exploiting the 
transformation properties in terms of spherical tensors, 
as will be shown next. 



7 



Sequences of Interaction Propagators 

Clearly, a palindromic sequence of propagators 

U s U s -i---U 2 U 1 U 1 U 2 ---U s - 1 U a (29) 

is locally invertible, if cither each component Uk or at 
least one partitioning of the sequence is locally invertible. 



Table III: Inversion of Spherical Tensors by Joint z-Rotations 
rotation angle <f> inverts interactions of quantum order ±p 



7T 

tv/2 
tt/3 



n/r 



1,3,5, ••• ,2q+l < j 
2, 6, 10, • • • , 4q + 2 < j 
3, 9, 15, • • • , 6q + 3 < j 



r, 3r, 5r, 



,r(2q + l) < j 



Multi-Qubit Interaction Hamiltonians of p-Quantum 

Order 



As usual in the treatment of angular momenta in spin-j 
representation (where j — n ■ j' may sum the spin quan- 
tum numbers of n identical, i.e. permutation symmetric 
single spins-j' to the group spin-j) one defines a rank-j 
spherical tensor Tj : , n of order m by the transformation 
properties under rotation by the Euler angles {a,/3,7} 



m'=-j 



(30) 



where the elements 



t->(7') i n \ i- ft — i— <jW5 —i^- <r {j] — i^er ' 



z J, m) 



(31) 



constitute the full Wigncr rotation matrix. An equivalent 
definition of the spherical tensors via the Pauli matrices 

or angular momentum operators {J x , Jy, Jz} = i su(2) 
in spin-j representation uses the commutation relations 

\Jx i iJyi Fj 7n j — [(/ , F^rn\ 



= v/j(j + l)-m(m±l) r iim±1 (32) 

and establishes the relation to the algebra sl(2, C ) rep- 
resented by {J+, J_,</ 2 }, as will be further illustrated 
below. 

Now specialise D^(a,f3,j) to D«' (0, 0, 0) by setting 
a = /3 = and y = <f>. Moreover, identify m = mf 
with the quantum order p of the interaction Hamiltonian 
H = Tj^p to obtain the eigenoperator equation 



(0,0,0) T iiP D«(0,0 



e" 1 ^ T iiP . (33) 



f G {x, y, z; +, — } by the symmetric sum over n qubits 
with in the £ th place written again as 

n 



£=1 



(34) 

Therefore, we may also envisage (0,0,0) as a /o- 
caZ z-rotation by an angle 4> acting jointly on all the n 
qubits, i.e., 0^(0,0,0) =: K((/),F Z ) G SU(2) <8 ". Thus 
one arrives at the two identical formulations 



(35) 



Clearly, inverting Tj tP to — Tj tP requires e~ lp ^ = — 1. By 
the Fourier duality between the quantum order p and the 
phase 4> one readily finds the results given in Table ITTT1 
A joint local z-rotation of angle n/r (with a fixed r = 
1,2,3,...) simultaneously inverts all the rank-j tensors 
of different quantum orders p given in the same row of 
the tabic. Thus also any linear combination of interaction 
tensors of quantum orders ±p = r{2q+\) can be inverted 
at a time for all q = 1, 2, 3, ... . 

Obviously, interaction Hamiltonians represented by 
spherical tensors X^o of order p = cannot be inverted 
by z-rotations K(cj), F z ), but may possibly be inverted by 
local unitaries operating on other than the z-axes. How- 
ever, rank-0 tensors T^o transforming like pscudoscalars 
such as e.g. the Heiscnberg XXX interaction Hamilto- 
nian, which in two spins is proportional to 



T, 



0,0 



2V3 



(a x ® a x + a y (g> a y + a z ® a z ) , 



(36) 



cannot be inverted at all. Not only does this hold for 
local unitaries, but for any similarity transform by some 
X G GL(2"), since the non-zero eigenvalues of To,o do 
not occur in pairs of opposite sign {vide supra). 



Inversion under Individual Local z-Rotations 



Inversion by Joint Local z-Rotations 

Note that the commutation relations for su(2) or 
sl(2, C) do not change, if one replaces the J„ where 



Since in a single spin qubit, Ad^ iZ has the eigenop- 
erators J„ G {1, J z , J+ , J_ } associated to the respec- 
tive eigenvalues e~ lp "^ G {1, 1, e~ p +^, e~ p -^}, it is easy 
to generalise the previous arguments to the case of 
z-rotations on n qubits — but with individually differing 
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rotation angles on each spin qubit <p±, <f>2, • ■ ■ , 4>i, ■ ■ • , 4>n- 
Now consider a Hamiltonian H taking the special form 
of a tensor product of Ad^, iZ eigenoperators on each spin 
qubit £ = 1 , . . . , n according to 



H 



J (1) 



4? 



■ ■ ® 49 ® ■ • • ® 4 n) 



(37) 



with independent i>i € {1, z,+,— } on each spin qubit. 
Then i? is clearly an eigenoperator to individual local 
z-rotations K(<j> x , ... ,<j> n , F z ) G SU(2)®" by virtue of 



Ad K 



(0i, ... ,0„,F Z ) 



_ g— »(Pi^i+ 



(38) 



So if is inverted if there is a set of rotation angles {<fie} 
with 



E 



±7r (mod 27r) 



(39) 



which is the case if there is at least one spin qubit I giv- 
ing rise to an interaction of quantum order pi = ±1. 
Moreover, a linear combination of such Hamiltonians 

m 

Hs '■= C \H\ is jointly invcrtible by an individual 

A=i 

local z-rotation K(4>\, . . . ,<fi n ,F z ) , if there is at least 
one consistent set of rotation angles {4>e} simultaneously 
satisfying for all the components H\ 



±7r (mod 2tt) 



which expresses the linear system 



/Pu 

P21 


Pl2 ■ 
P22 ■ 


■ Pln\ 

■ P2n 


(<t>A 

4>2 




P31 


PZ2 ■ 


■ PZn 






\Pml 


Pm2 ■ 


Pmn/ 


\<t>n) 





±7T (mod 2ir) 
±7r (mod 2ir) 

\±7r (mod 2tt)/ 



(40) 

Note the signs on the RHS may be chosen independently 
in 2 m ways with every choice forming a system of 
m linear equations in n variables. Therefore, if the 
vector of any of the combinations 7r (±1, ±1, ••• ,±1)' 
can be expanded in terms of the column vectors of 
P := (px,e) with real coefficients, then Hy, is locally 
invertible by z-rotations. In the special case of m = n 
and P non-singular, there always is a consistent set 
of individual rotation angles for inverting for any 
choice of signs. For simplicity, we will drop the index 
£ in Hy: henceforth writing H for Hamiltonians locally 
invertible by individual ^-rotations. 



Corollary 1 Let H be locally invertible by individual 
z -rotations on each qubit. Then the following hold. 

(1) Any Hamiltonian H on the local unitary orbit 
AdK (H) of any such H generates a one-parameter 
unitary group that is locally invertible of type-I. 



(2) In turn, any type-I locally invertible Hamiltonian 
H is on a local unitary orbit of some H . 

Proof: (1) is obvious. (2) follows since every local uni- 
tary K is locally unitarily similar to a local z-rotation K z : 
K = KK Z K~} => KHK- 1 = -H <s> K z HK~ l = -H 
where H := K~ X EK. ■ 
Thus local invertibility by z-rotations can be looked upon 
as the normal form of the problem. 

Moreover, in order to see the link to the root space 
decomposition of the semisimple Lie algebra sl(2, C ) in 
spin-j representation, take the derivative of Eqn. 1351 



d 



6=0 Ac W, 



J:l> 



d_ 



1=0 



-iad kz (T j>p ) 



-i P T. 



3,P 



WTj- (41) 



(42) 



(NB: the minus sign on the left side of the last identity 



is due to the convention K((j>, z) 



D —i<fik z 



UH 



analogous to 



imposed by Schrodinger's equation.) 



U :-- 

Note the tensors are eigenoperators to ad^, of the joint 
local z-rotations as anticipated in Eqn. 1321 



Algebra II: Root-Space Decomposition 

To fix notations, let q be a complex semisimple Lie 
algebra and let go be a Cartan subalgebra of g, i.e. a 
maximally abelian subalgebra such that for all H E go, 
|34| the commutator superoperators ad# are simultane- 
ously diagonalisablc. Then the root-space decomposition 
of g with respect to go takes the form 



= 0o©0i 

ctjtO 



(43) 



where the root-spaces g Q (with a ^ 0) are the non-trivial 
simultaneous eigenspaccs 



9a ■= {g G 9 I ad H (.g) = a(H) g} 



(44) 



The corresponding non-trivial a's are called roots of the 
decomposition. They are elements of the dual space g^ 
of linear functionals on go- 

In the following, we consider the complex semisimple 
Lie algebra sl(N, C ) as the complcxification of the real 
Lie algebra su(N). Define Eij as a square matrix differ- 
ing from the zero matrix by just one element — the unity 
in the j th column of the i th row. Moreover, let go be 
the set of all diagonal matrices in sl(N, C ) and define 
ei(H) := Ha for all H € go- Then for every H G g , the 
Eij are simultaneous eigenoperators of the commutator 
superoperators ad# with eigenvalues depending linearly 
on H 



win (Eij) — (Ha — Hjj) E^ 

= ( ei (H) - ej (H)) E k 



(45) 
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Thus the root-space decomposition of sl(N, C) may be 
rewritten as 



where the box sizes coincide with 1. 



= So © C E i:j 



(46) 



Furthermore, if g' is a compact real semisimple Lie al- 
gebra with maximally torus algebra t', then the complex- 
ification t of t' gives a Cartan subalgebra of the complex- 
ification g of g'. For example in the case of su(iV), 



t':={i diag(0i,0 2 



X> = o> 



may be chosen as maximal torus algebra. Then 
t := t' + it' 



(47) 



(48) 



i.e. the set of all complex diagonal matrices forms the 
Cartan subalgebra of sl(N, C ). 

Now Ean. l45l shows that an Eij with i ^ j can be sign- 
inverted provided (Ha — Hjj) ^ 0. The generic case, 
joint and individual local z-rotations are specified next. 

Proposition 1 In a system of n spins-^, for the single- 
element matrices Eij with i ^ j the following hold: 

(1) to any E^ there is an element of the Weyl torus 
taking the form T = exp(— i diag (6i, 6*2, ■ ■ ■ , #jv)) 
so that AdT(Eij) = —Eij; 

(2) any matrix Eij can also be sign- inverted by a single 
local z-rotation; 

(3) in contrast, by a joint local z-rotation on all the 
n spins-^, an Eij can only be sign-inverted if for 
its indices i,j the reductions by 1 written as binary 
numbers (i — 1)2 and (j — 1)2 do not have the same 
number of 0's and l's (irrespective of the order). 

Proof: First, note that although s[(2",C) comprises 
the generators of all special unitary propagators, its 
maximally abelian subalgebra go can — without loss of 
generality — always be chosen such that it includes the 
generators of all the local z-rotations. They suffice to be 
considered, since the Eij are simultaneous eigenopcra- 
tors to all elements in go. 

(1) Obviously elements in the Weyl torus algebra t can 
be chosen such that 9i — 6j 7^ 0. 

(2) Clearly any off-diagonal element E^ in the boxes of 
1 □' 

P -1 

zero root (e^ — ej){H) = Ha — 
for H <E> 1 and 1 <g> H. Likewise any off-diagonal element 
Eij is in one of the boxes of the following embedded Pauli 
z-matriccs 



the block matrix H := 



is associated with a non- 



Hjj . The same holds true 



0~£z 



i(£-l) 



1 □ 
□ -1 



1 



®(n-£) 



(set If 



1). 



Let the index run I = 1,2, ... n to see that in fact 
every off-diagonal element E^ can be associated with 
some i, which implies any E^ can be sign-inverted by 
at least one local z-rotation on some single spin qubit 
£. (Due to permutation symmetry, for £ < n there 
are off-diagonal Eij with non-zero roots even outside 
the boxes; they add further options of choosing a qubit I.) 

(3) For Eij, let (i - 1) 2 =: £fe=o 2 k b k and (j - 1) 2 =: 
Sfc=o 2 fc &fc define the ri-digit binary representations of 
the indices reduced by 1. If (i — 1)2 and (j — 1)2 have the 
same number of 0's and l's, we will show that E^ belongs 
to a zero root of F z , i.e. adi?_(i?ij) = 0. In Appcndix-B 
we gave a general formula for the matrix elements (F z )a. 
So 

n—1 n—1 

(F*)u ~ (F z )jj = i ]T(-l) bfc - \ £(-l) 6 '* (49) 



fc=0 



fe=0 



vanishes, if and only if an equal number of terms (— 1) 
and (— l) 1 appears in both n-term sums as claimed. I 

Corollary 2 As the maximally abelian algebra Qq of both 
s[(2,C) and su(2) in the representation of n spins-j' can 
always be chosen such as to comprise the generators k z 
of K((j), F z ) G SU(2)®" bringing about local z-rotations 
jointly on all spins, the tensors of order p are associated 
to the root space elements Eij o/s[(2,C) showing the 
eigenvalue (e^ — ej){F z ) = p. 

Allowing individual z-rotations on each qubit again, 
one finds 



Ad 



(Eij) 



-i(pi<f>i+ ■■■ +p„ 



E, 



(50) 



because any E^ can be written as a tensor product of the 
single-element two by two matrices {J a , J 13 , J + , J_} := 
{(58), (§?), (85), (?8) } associated with the eigen- 



values 
and p± = ±1. 



for v G {a, (3, +, — }, where p a 



VP 



(0 



Examples: 

(1) The single-element Weyl matrix -Eos, 15 belongs to the 
zero-root (e, — ej)(F™) = p = since for all n > 4 the 
binary representations (according to Proposition 1.3) end 
with 72 = 0111 and 142 = 1110 having the same number 
of 0's and l's. Thus it cannot be sign- inverted by joint 
local z-rotations, whereas by being off-diagonal it can 
always be sign-inverted by an individual local z-rotation. 

(2) In contrast, for -^47,11 one finds by Ean. 1491 and the 
binaries 46 2 = 101110 and 10 2 = 001010 that p = (-4)|. 
So it can be sign-inverted by a joint local z-rotation with 
rotation angle <fi = ? in accordance with Tab. IIIII 

Given the relation to the transformation properties of 
spherical tensors, it is easy to analyse type-I local invert- 
ibility of linear combinations of root space elements E^ 
under joint or individual z-rotations. 
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p = p = ±1 p = ±2 p=±3 p=±4 

Figure 3: (Colour online) Matrices Eij constituting the rank-4 spherical tensors Tn, v of order p (here in a system of 4 spins-^, 
p 6 [—4, +4]). The non-zero elements in the Eij are marked ■ for p > and ■ for p < 0. The T4, p are eigenoperators to the 
Weyl torus element F z 6 t according to ad^ (T4,p) = P ?4,p with the eigenvalues being the quantum orders p. The constituents 
of the tensors are the single-element Weyl matrices Eij sharing the same eigenvalues ad F z (Eij) = p Eij. According to Tab, IIIl| 
a Hamiltonian comprising elements that transform like Ta,± p can locally be sign-inverted e.g. by a joint ^-rotation of angle 
n/p. 



Proposition 2 In a system of n qubits, a linear combi- 
nation of single- element matrices 



(A) 
3 



(51) 



with i j and c\ G C is sign-invertible by an individual 
local z-rotation K(cj>i, ... ,<fi n ,F z ) , if there is at least 
one consistent set of rotation angles {4>i\£ = 1,2, ... ,n} 

simultaneously satisfying for all its constituents E^ 

n 

J~]p\,t • <j>t=±n (mod 2ir) , 
which coincides with the linear system in Eqn. \40\ 



Relation to Time Reversal and Cartan 
Decompositions 

As will be shown, the detailed discussion of the root- 
space decomposition in the previous section was in fact 
needed, and a mere Cartan decomposition does not de- 
cide about type-I local invertibility. 

Let g be a real compact semisimple Lie algebra and let 
the mapping 9 : g — ► g be any involutive (Lie algebra) 
automorphism. Then 9 defines a Cartan-like decomposi- 
tion g = ig © pg of g, where tg and pg are the respective 
+1 and —1 eigenspaces of 9, i.e. 

6{X) = X for all X e l e (52) 
6{X) = -X for all X e p e (53) 

ensuring the standard commutation relations 

[loM C tg (54) 

[te,Pe] Q Pe (55) 
[pe,Pe] C tg . (56) 

In su(2"), one may choose the so-called concurrence Car- 
tan involution |27j 



(X) := (-ia y r n X*((-ia y r n y 



where 9qc takes the form of the bit flip operator and 
thus relates to time reversal. Bullock et al. classi- 
fied Hamiltonians iH € pec as symmetric with respect to 
time-reversal and those in ice as anti-symmetric. Since 



{-i<Ty) 



-inF v 



this representation of the Cartan in- 



volution is equivalent to a local y-rotation acting jointly 
on n qubits following complex conjugation. Due to the 
latter, the Cartan involution 9 is unphysical (as also 
pointed out in ref. (2^) and it is thus distinct from the 
local unitary operations discussed here. 

Note that in su(4) ice coincides with the algebra 6 
generating the local unitaries K := SU(2)® 2 , which is the 
reason for our notation, whereas in su(2 n ) with n > 3 this 
is no longer true, since 6cc comprises m-linear interaction 
Hamiltonians with m odd, while pec encompasses those 
with m even. 

As described above for the simple case of two qubits, 
the pair interactions (m = 2) Hzz,Hxy £ Pec are lo- 
cally type-I invertible, while ifxxx, Hxyz S Pcc are not. 
Yet in two qubits tec = t (s.a.), so all the elements in 
?cc are — by definition — type-I invertible, while for 3 and 
more qubits, tec generically contains type-I invertible in- 
teractions (e.g. zzz) as well as non-invertible ones (e.g. 

-^xxx-f 2yyy+3zzz ) ■ 

Hence, a Cartan-type decomposition into time-reversal 
symmetric and antisymmetric subspaces does not decide 
whether an interaction is locally invertible or not. In 
SU(A r ) with TV := 2™, also for other standard choices of 
the Cartan involution, such as [2^ 



ui(X) 

MX) 



x* 

JnX*J n 

Ip.qX Ip.q 



with p = q = N/2 in the definitions 



Jn 



(57) 



l P,Q 



1 

ljV/2 

, o 

-X 



N/2 



(58) 
(59) 
(60) 



(61) 
(62) 



11 



Table IV: Relation to Cartan Decomposition for Different 
Choices of Involution 



-^interaction Type-I 6qC &A I &A II III 

(number of qubits m) invertible 



r^auli 


V 1 

A 1 


+ 


t 


P 


t 


p 


matrices 


A r 1 


+ 


<> 
t 


fc 


t 


p 




Zl 


+ 


I 


P 


P 


t 


m = 2 


Hzz 


+ 


P 


P 


P 


t 




a xx 


+ 


P 


P 


t 


p 




Hxy 


+ 


P 


P 


t 


p 




Hxxx 




P 


P 


tup 


tup 




Hxxy 




P 


P 


tup 


tup 




Hxyz 




P 


P 


e u p 


e u p 


m = 3 


a zzz 


+ 


f 


P 


t 


t 




a xxx±yyy 


+ 


t 


6 U p 


t 


p 




H X xx-\-yyy-\- zzz 


+ 


t 


t U p 


t 


tup 



Note: 

8cc and 9 a ii are equivalent up to non-local permutation; 
the same holds for 9 a ii and 8a hi in the case of p — q. 



the decomposition into ie and pg does never completely 
agree with the subdivision into locally invertible and 
non-invertible interacion Hamiltonians iH lnt £ su(2") as 
shown in Tab. IIVI Rather in the general case, it takes 
the more specific patterns derived from the root space 
decomposition as described in the previous section. 

Gradient Flows for Type-I Inversion by Local 
Unitaries 

Finally there is — fortunately — a convenient numerical 
solution to the decision problem whether a given Hamil- 
tonian H generates a locally invertible unitary, which is 
particularly helpful in cases where algebraic assessment 
is tedious. It recasts the problem to the question whether 
the minimum of the distance 

A(K) := \\KHK- 1 + H\\ 2 = || Ad K (H) + H\\ 2 (63) 

over all local unitaries K £ SU(2)®" is zero or not: 
clearly, the norm ensures that A (if) = if and only 
if AdK(-ff) = —H, which means [35| 

f(K) := Retr{K H K^ 1 (- H)} (64) 

shall attain a global maximum that has to coincide with 
the upper bound for hermitian H reaching equality in 
fmax(K) < \\H\\2- Whether this limit can be reached 
by local unitaries may rcadlily be checked numerically. 
To this end, one may devise a gradient flow along the 
lines of Ref. [2^, [3(| where, however, the gradient in the 



1 - 




50 100 150 

iteration 



Figure 4: (Colour online) Gradient-flow driven local inversion 
of different Heisenberg interaction Hamiltonians: (a) the ZZ 
interaction on a cyclic four-qubit topology Ci can be locally 
inverted, (b) the ZZ interaction on a cyclic three-qubit topol- 
ogy C3 cannot be inverted locally, (c) nor the XXX isotropic 
interaction between two qubits. 



1 - 




25 50 75 

iteration 



Figure 5: (Colour online) Gradient-flow driven local inver- 
sion of Heisenberg interactions of different quantum order on 
a cyclic three-qubit topology C3: (a) the double-quantum in- 
teraction X(— X) can be locally inverted, while (b) the zero- 
quantum analogue XX cannot. 

tangent space has to be restricted by projecting it onto 
the algebra of local unitaries t generating SU(2)® n . As 
will be described elsewhere, establishing convergence and 
appropriate step sizes of the iterative numerical scheme 
can be handled on a very general level. 

In the present context notice the gradient flow on the 
local unitaries takes the form 

K = gradf(K) = PtdiKHK- 1 ),!!}) K 

(65) 

= -P E (ad H o Ad K (H))K, 
where P{ denotes the projection onto the subalgebra £ 
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of generators of local unitarics K = SU(2)® n . The flow 
clearly reaches a critical point if already the entire com- 
mutator vanishes 

[Ad K (H),H}=0 . (66) 

For hermitian if, this is the case, for instance whenever 

Ad K (H) = e ±lpTr H (with p = 0,1, 2,...) , (67) 

which means H is an eigenoperator to Adx(-ff), i- e -; 
vecH is an eigenvector of [K* ® K). Eigenvectors H 
to the eigenvalue +1 lead to global minima of f(K), 
while global maxima are reached by eigenvectors H to 
the eigenvalue — 1. 

In Figs. El and we give some examples. Let H be 
normalised to ||-ff|| 2 = 1. If ^{KHK^H} = -1 can 
be reached, the interaction Hamiltonian is locally invert- 
ible as in the case of the Heisenberg ZZ interaction in 
a cyclic four-qubit coupling topology (which clearly is a 
bipartite graph) , while in the cyclic three-qubit topology 
(obviously not forming a bipartite coupling graph) or in 
the case of the isotropic XXX interaction it is not. 

Relation to Local C-Numerical Ranges 

The C-numerical range is well-known |3lT ] to consist of 
the following set of points in the complex plane 

W(C, A) := {tr (C^UAU- 1 ) \ U e SU{2 n )} . (68) 

In Rcf. [32, we defined as local C-numerical range its 
subset 

Wi oc {C, A) := {tr (C f KAK' 1 ) \ K e SU{2f n } . (69) 

In view of locally reversible Hamiltonians, things spe- 
cialise to C = —H = —A. Normalising again to \\H\\ = 1, 
a locally reversible Hamiltonian H clearly requires — 1 € 
Wi oc (— H, H). This has just been exemplified by the nu- 
merical examples in the previous section. Moreover, be- 
ing a linear map of the local unitary orbit, the local C- 
numerical range is connected. For locally reversible H, 
one finds the real line segment [— 1; +1] = W\ oc {— H, H), 
whereas in Hamiltonians that fail to be locally reversible, 
the line segment falls short of extending from +1 (which 
trivially always can be attained) to —1. 

With these observations, the different aspects may be 
summed up. 

Synopsis on Type-I Inversion 

Corollary 3 (Local Time Reversal) For an interac- 
tion Hamiltonian H = with ||-ff|| 2 = 1 the following 
are equivalent: 

1. H is locally sign-reversible of type- 1; 




Figure 6: Venn diagram showing that simple criteria of non- 
zero eigenvalues in pairs of opposite sign, adl(H) = H, and 
Ijj n t 7^ {} fail to decide type-I invertibility as explained in 
the text. The numbers in the sets refer to the examples listed 
in Tab. El 

2. its local C-numerical range comprises — 1: 
-16 Wi oc (-H,H); 

3. its local C-numerical range is the real line segment 
from -1 to +1: W\ oc {-H,H) = [-1; +1]; 

4. 3K 6 SJ7(2)®" : 

\\KHK- 1 + H\\\ = ^ Ad K (H) = -H 

5. H is locally unitarily similar to a H with 
Ad Kz (H ) = -H ; 

6. let Q = Qo © ® C Eij be the root-space decompo- 

sition of sl(N, C); H is locally unitarily similar to 
a linear combination of root-space elements to non- 
zero roots 

A=l 

satisfying a system of linear equations 
y^py • 4>t = ±7r(mod 2ir) 

in the sense of Eqn. \4t\ where the p\,e can be in- 
terpreted as the quantum orders of the constituting 
spherical tensor elements. 

Proof: The equivalence of (1) with statements (2) 
through (6) was of course already proven in the respec- 
tive sections. 

Moreover, one finds (1) (2): obvious; (2) (3): con- 
nectedness of W\ oc {C,A); (3) =>• (4): obvious; (4) =*> (5): 
Corollary 1; (5) =>■ (6): Corollary 1, Proposition 1 and 2; 
as well as Corollary 2 for the interpretation as quantum 
orders; (6) (1): Proposition 2. ■ 
However, the simple necessary criteria of (i) non-zero 
eigenvalues occuring in pairs of opposite sign, (ii) the in- 
tersection of the orthocomplcmcnt to the invariant sub- 
space with the generators of local unitaries not being 
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Table V: Examples of pair interaction and multi-qubit in- 
teraction Hamiltonians used to show the coverage by simple 
type-I inversion criteria in Fig. |S| 



Example 


Hamiltonian 


1 


zz 


2 


xx + yy 


3 


xxl — yyl + xlx 




-yly + lxx - lyy 


4 


xxll — yyll — lxxx + lyxy + lyyx + lxyy 




—xlxx + ylxy + ylyx + xlyy 


5 


xxll — yyll + xlll — lxxx + lyxy + lyyx + lxyy 




—xlxx + ylxy + ylyx + xlyy 


6 


zll — xxx + xyy + yxy + yyx 


7 


xxl + yyl + zzl 




— (lxx + lyy + lzz) 


8 


xxl + yyl + xlx 




+yly + lxx + lyy 


9 


zzl + zlz + lzz 


10 


zz + zl + lx 


11 


zz + zl + lz 



empty ijj fit 7^ {}, as well as the sufficient condition (iii) 
of the double commutator reproducing the Hamiltonian 
in question, ad^(-ff) = H, fall short of giving a conclusive 
decision on type-I invertibility, sec Fig. 



II. POINTWISE LOCALLY INVERTIBLE 
PROPAGATORS 

Propagators that are not jointly invertiblc by a local 
unitary for all t € R (together with the entire one- 
parameter group generated by their Hamiltonian) may 
still be pointwise locally invertible at certain times r. So 
the task in this section is the following: given some t > 0, 
determine whether there is a pair \K\,K 2 ^ K% } C 
K := SU(2)®" so that 

K x e'" H K 2 = e +lTH 
(K\ <X> K x ) vec{e~" H ) = vec{e +lTH ) . 

Remark 2 Note that type-II invertibility only arises 
upon restriction to local operations K\,K<i G K, because 
to any Uq G SU(iV) there is a trivial pair U\,U2 G 
SU(iV) with U 1 U U 2 = U^ 1 (e.g. U t = U^ 2 ,U 2 = 1), 
whereas with K\ , G K there is no such trivial generic 
solution unless Uq G K. 

Corollary 4 Let H generate a one-parameter unitary 
group U := {e~ itH \t G R , iH G su(iV,C)} that is lo- 
cally invertible of type-I. Then 



1. the generic elements of the left and right cosets ¥3A 
and WK are type-II locally invertible; 



2. in turn, every Hamiltonian that is type-II invert- 
ible is an element of a coset YJA or UK., where U 
is some one-parameter unitary group that itself is 
type-I invertible. 



Therefore type-II invertiblc propagators are a natural 
extension of the typc-I invertiblc unitary one-parameter 
groups. In turn, however, the decision problem whether 
a given propagator is type-II invertible is generally quite 
complicated so that we will devise a coupled gradient flow 
on two local unitaries for solving it numerically. Yet a 
number of cases can be treated algebraically by analysing 
the symmetries of the matrix representation of the uni- 
tary propagator to be inverted. 

Since these symmetry considerations extend beyond 
the representation of unitary matrices, we will ask 
whether an arbitrary given matrix can be mapped to 
its hermitian adjoint by a superoperator of the form 
(K\ <g) K\) with local unitary K 1 ,K 2 (cp. Eqn.[TO|). To 
this end, one has to maximise the coincidences between 
{K\ (g) K x ) and the adjoining superoperator denoted Adj 
that takes its argument to the hermitian adjoint (i.e. the 
complex conjugate transpose). Clearly, there is no local 
unitary {K\ ® K x ) that fully matches with Adj as this 
would be a universal inverting operator. However, there 
are classes of partial overlaps, where the lack of coinci- 
dence enforces a symmetry in the matrices to be adjoined. 
These will be analysed in detail in the following. 

Because Adj has no matrix representation over the 
field of complex numbers, we turn to the real domain. 
With Mr and Mi m denoting the respective real and 
imaginary parts of an arbitrary complex matrix M , one 
obtains a convenient representation of M as a real vector 
by virtue of the faithful mapping 

M i ► vcc(M Re ) © vcc(M Im ) . (71) 

[Note that this representation shows less redundance 
than the usual M h-> ( ^ ) ,] 

In this notation, the adjoining superoperator docs have 
a real matrix representation defined via 

Adj H (vec(M Re ) © vec(Mi m )) = vec(M^ c ) © vec(-M I * m ) 

(72) 

such as to take the form 

Adj R = (J_J) , (73) 
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by virtue of the transposition supcroperator T, which e.g. 
for the above representation of a M £ Mat4(C ) reads 



z 1 



T : = 



\ 



\ 



(74) 



1/ 



Likewise, for the local unitary transform 
KxMK-2 = Kvec(M) with K := K\ © K x one gets the 
corresponding real representation of the superoperator 
via 



Km vec(M) R := 



-Ki 
Ki 



vec(M Re ) 



vec(Mi m )) . 
_ (75) 

Comparing the structure of Km here and Adj R (in 
Eqn. 173(1 immediately shows that for maximal coinci- 
dence the imaginary block K\ m within Km has to vanish, 
because the row and column norms are limited to unity 
in (local) unitaries. One may readily express the utmost 
possible overlaps of Km and Adj R by taking the clement- 
wise Hadamard product as the coincidence matrix C 



Cp 



Km © Adj ! 



ATrc -K lm 
Kim Kr, 



C 
-C 



© 



( f I) 
\ — T / 



(76) 



where C reads, e.g. in the case n = 4 

/+a 



C:=± 



+d 



= d 



+d 



+b 



V 



= b 



■ +a J 

(77) 

and in which either a or b or c or d is unity. Thus 
in the case n = 4 one finds four subtypes of max- 
imal overlap, termed A,B,C,D henceforth. Accord- 
ing to the possible choices of signs, each of them oc- 
curs in four sign patterns expressed by the indices v £ 



++++! - 



A+,,A,+_,A. 



_} and analogously for 
subtypes B,C,D. 

For instance, let v = A-\ |_, then the local unitary 

Km for maximal overlap with Adj R shows the following 



Table VI: Subtypes and Sign Patterns in Type-II Inversion 



# spin qubits # subtypes 



# sign patterns 



non-zero block K^ e : 
/ + 



K 



(A+ — + ) 
Ro 



(78) 



7 



where the elements ± stand for ±1, while * are unavoid- 
able non-zero elements enforced by Km being a local uni- 
tary. They do not contribute to Adj R , in contrary, they 
bring about unwanted actions on the argument, i.e. the 
matrix M £ Mat „ (C ) . Together with the lacking ele- 
ments for full overlap with T, they require the following 
symmetry in the matrix argument 



Mi 



MR 








Ml4 















±M 2 * 3 







t«; 4 








M Ra 

M 44 



(79) 



in order to fulfill Km vec(M) = vcc(M^) as desired. 

For the sake of completeness in the case of n = 4, 
we give the remainder of constituents in the subtypes 
A,B,C,D as well as the associated sign patterns in 
Appendix-C and D. 

The structures of the pertinent block matrices T, Adj K 
and hence C are easily scalable to larger n: in Tab. [23 
we give the number of subtypes of coincidence as well 
as the number of different symmetry subtypes and sign 
patterns in the matrix arguments M £ Mat„(C) with 
growing number of dimensions. 



Type-II Inversion via Coupled Gradient Flows on 
Two Local Unitaries 



In the general case, one may conveniently restate the 
problem of pointwisc local invertibility to the question, 
whether for a fixed non-zero r £ K there is a pair 
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20 30 
iteration 



Figure 7: (Colour online) Gradient-flow driven local inversion 
of U(t) := exp{-if H} with H = (ff z ®l+l®a z +a z ®a z )/2 

(a) by coupled gradient flows on independent Ki and K2 and 

(b) by a gradient flow with K 2 = K± . 



K 1 ,K 2 GK = SU(2)® n so that 



K ie - iTH K 2 = e +irH 



^\\K 1 e- lTH K 2 -e +irH \\ 2 = 



(80) 



Then one may devise a coupled gradient flow on two local 
unitaries simultaneously in order to minimise 



f(K 1 ,K 2 )=Rc tT{K ie -" H K 2 (-e- lTH )} 

by (writing U := eT %TH for short) 

K x = grad/(X 1 ) = P t (KxUK^-U)) K x 
K 2 = gradf(K 2 ) = P t (K^-U^U) K 2 



(81) 

(82) 
(83) 



Again, if ^ Re tr{K ie - tTH K 2 (- 



,-irH 



)} = —1 can be 
reached, then U(t) = e~ %TH is locally invertible at the 
point t. Examples are shown in Fig. [7| 



closed algebraic form. Hamiltonians generating type-I 
invertible evolutions are locally unitarily similar to those 
invertible solely by local z-rotations, which thus can 
be regarded as normal form. Taking the differential, 
we showed their components E^ relate to the non-zero 
roots of the root-space decomposition of sl(N, C ) via 
adfe_ (Eij) = (e.i — ej)(Eij), where k z are the generators of 
local ^-rotations. Moreover, in the special case of joint lo- 
cal z-rotations generated by F Zl the non-zero roots were 
further shown to relate to the spherical tensors T J;P of 
non-zero quantum order p. For Ising ZZ-coupling inter- 
actions as well as for Heisenberg XY interactions to be 
locally invertible of type-I, their coupling topology has to 
take the form of a bipartite graph. An exception is the 
Heisenberg X(— X) interaction, which is typc-I invert- 
ible on any coupling topology, while Heisenberg XXX, 
XXY, and XY Z interactions are not locally invertible 
at all because they relate to rank-0 tensors, and their 
non-zero eigenvalues do not occur in pairs of opposite 
sign. 

The pointwise invertible quantum evolutions of type- 
II are generalisations of those of typc-I. They consist of 
cosct elements K.U and UK., where U are type-I invertible 
one-parameter groups. Here we have caracterised type-II 
invertible propagators by the symmetries of their matrix 
representations . 

Finally, in view of convenience in practical applica- 
tions, we have devised gradient-flow based numerical 
checks to decide whether 



mm 

K 1 ,K 2 GSU{2)i 



\K 1 e- ltH K 2 -e +UH \\l = Q , 



i.e. whether a propagator is locally invertible of type-I 
or type-II or not. 



III. APPENDIX 



A. Classification 



Conclusion 

Generalising in the sense of Halm's spin echo, we have 
characterised all those effective multi-qubit quantum in- 
teractions allowing for time reversal by manipulations 
confined to local unitary operations. The evolutions 
generated by these interaction Hamiltonians can be re- 
versed and refocussed solely by local unitaries. To this 
end, we have given a number of necessary and sufficient 
conditions in terms of geometry, eigenoperators, graphs 
of coupling topology, tensor analysis and root-space de- 
composition. Moreover, we have classified locally in- 
vertible evolutions into two types. Type-I consists of 
one-parameter groups generated by Hamiltonians H that 
are eigenoperators of local unitary conjugation associated 
with the eigenvalue —1, i.e. A&k{H) = —H. We have 
shown how to construct the corresponding eigenspace in 



First, we prove Lemma from the introduction: 
Either e~ ltH is 

1. not invertible by local unitaries at all, or 

2. it is trivial and self-inverse, or 

3. it is type-I invertible in the sense 3K e SU(2)® n : 
KHK'- 1 = -H so Ke~~ ltH K^ 1 = e +itH jointly for 
all t G M , or 

4. it is type-II invertible such that at some (but not 
all) points r in time Kie~ lrH K 2 = e +lrH with 
K u K 2 e SU(2)® n and K 2 ^ K^ 1 . 

Proof: By the series expansion of the exponential one 
finds the obvious equivalence 

KHK- 1 = -H 

& Vf G R : Ke~ ltH K- 1 = e +itH . 
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Its logical negation 

KHK~ l ^ -H 
^> 3t e 



Ke 



itH 



K- 1 e 



(85) 



comprises the following trivial cases 

1. Vt S K : Ke~ itH K^ 1 ^ e +ltH , so either e~ UH is 
not locally invertible at all, or 

2. 3t : Ke~ irH K^ 1 = e +lrH , while for all other 
t =/= t (with exceptions of measure zero due to pe- 
riodicity) Ke~ itH K~ x ^ e +ltH while (KHR- 1 ) ^ 



H. This can only hold, if K = 1 and e 



rH 



is 



self-inverse. 

Otherwise, if the affirmative (Ean l84|) is true one has 



Proof (induction): 

For n = 1 one has: (i - 1) 2 = 2°6 G {0,1} so 
= iEfro 1 ^-!) 6 * S ivin S (^)n = I = -(^) 22 - 

In order to proceed from n — > n + 1 we show that 
with (-F 1 ™)^ being given one finds for the new index 
i' :=2"C + !€ {1,2,3,..., 2™ +1 } 

(^ n+1 W =(l®F 2 ") i v + i(-l) 6 " • (88) 

Use 

F 2 " +1 = 1 2 $ F z ™ + I, ® l 2 n 

= diag ((i?) n , . . . , (F«) u , . . . , (i?W) 

+ ±diag(l,l,...,l;-l, -!,...,-!) 



3. Vt G R : Ke~ itH K^ 1 = e 



itH 



(89) 



Finally, we have to show that type-I and II are distinct to see that the last term adds | for il g {1, . . . , 2™} 



4. K ie - tTH K 2 = e +tTH with A" x ^ JCT 1 may hold 
pointwise for certain r, but not for all r G R . As- 
sume the contrary: Vt G R Kie~ ltH K 2 = e +itH 
with K 1 7^ ^2 an d define as commuting ele- 
ments of a one-parameter group U\ := e~ ltlH and 
U 2 ■= e~ lt2H to give U 12 := e - i( - tl+t ^ H . Then one 
has 

K 1 U 12 K 2 = U^ 2 
K X U X U 2 K 2 = U^Ui 1 = Ui 1 U 2 l 
{K\ Ui K 2 ) (K 2 1 U 2 K 2 ) = {KxUiK 2 ){K^J2K 2 ) 

K 2 l = K x , 

where the latter contradicts the assumption. 

These four instances prove Tab. [I] ■ 



and -\ for %' G {2 n +l, . . . , 2" +1 } = 2 n +i, in coincidence 
with b n taking the value or 1. ■ 



B. Explicit General Representation of F z 

Recall that the generator of a joint z-rotation on all 
the n spin-i qubits is defined as the diagonal matrix 

n 

F * ■= 2 E 1 (1) ®1 (2) ®' ■ ■®1 ( ^ 1) ®<tW®1^+ 1 )®- • 

(86) 

summing over the Pauli matrix cr^ on all qubits. When- 
ever it is necessary to express the total number of qubits, 
we write F™. Here we prove an explicit formula giving 
its i th diagonal clement for general n. 

Lemma 8 For (F?) u with the index i G {1, 2,3,..., 2"} 
calculate the n-digit binary representation for the reduc- 
tion by 1 as (i — 1) 2 =: Y^k=a 2 fc &fc- Then the i th diagonal 
element reads 

n-l 

ra« = IE(- 1 ) 6fc ■ ( 8? ) 

k=0 



17 



C. Subtypes of Pointwise Invertible Local Unitaries Subtype C has the forms 



For the case of two qubits, we give local unitary su- 
pcropcrators of different type of partial overlap with the 
adjoining superoperator 



K m vec(M)j 



K Ro -K lm 
Kim i?Ro 



>ec(M Ro ) © vcc(M lm )) . 

(90) 

In subtype A, the blockmatrix Kn, e within the above su- 
permatrix K m may take four different forms according to 
the indices 



/ + 



-4++++ = -A 



/ + 



A+-+ = -A. 



Subtype B comprises the forms 



/ 



\ 



/: 



C_| h — -C — h+ _ — 



V 

Subtype D includes the forms 



= -D. 



V 



B 4 



-B 



V* 



B 4 



-+ = -B-+ 



\ 
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D. Symmetries in the Argument 



According to the classes of partial overlap with the ad- 
joining supcropcrator, here we give the according sym- 
metries for the matrices M € Mat 4(C ) to be mapped to 
their adjoints by the corresponding local unitaries. 



Ma, 




C+ + + + 



Mn 


Mia 


MR 


M14 


M 2 i 


M 22 





M Rc 


mr 


M 32 











MR 












-l a A/* 

-i'at* 

A^l*2 



14 







Mii 



22 ■ 















1"M* 2 


n 










m; 3 








M* ti 


-l b M 2 * 4 


-1°M* 4 






M c+ . 



A/12 


l.f Im 
M 13 











MR 








±M 2 *i 


71 rim 
m 13 


ta/T 2 






M A+ _ 



MR 








A/14 





M 22 " 


A/23 








±M* 3 


MR 





TA/i* 4 








A/ Ro 
JW44 



M c _ 






Mi 2 


Af r<c 
lvl 13 





A/21 








mR 


< 








±M 2 *! 





MR 


TA/i* 2 














A/14 





M Re 


M23 





(1 


±A/ 2 * 3 


mR 





TM* 4 








A/ Im 

Ju 44 



Mr 



A/11 M12 Mj™ M14 
A/21 A/ 2 2 MR 
Ml: A/32 
MR 



M A _ 



M« 






Mi; 

Mil 





A/13 Mi. 

M 23 A/ 2 

MR A/ 3 

Ml! 



-(-1 Q ' 6 ) 






-l a A/* 
-l b A/i 

A^l*2 



14 



11 






A'2*i 
-l b M,*, 



-(-l a ' b ) 















l"A/* 2 


n 










M 2*3 








M^ 


-l b M* 4 


-1°M* 4 






M 



B+ + + + 



M f 



M f 



M E 



A/11 


A/ Rc 


A/13 


A/14 


MR 


(1 


A/23 


M 24 








A/33 


Jw 34 








A/ Ro 
M 43 






+(_!»•») 





-1°M* 

A/ 24 -l b A/- 

-l b M 23 M* 



11 



'1:1 




» flm 

M 12 


A/13 











M 24 








Af Ro 
Jvl 3i 


T A/* 3 


Jw 43 








f ° 


M Ro 

Jw 12 


A/13 




MR 










±A/*4 










V 


TM* 3 


M 43 






Mm 



M 



M 



D- 



Mu 


A/12 


A/13 


A/ Rc 

Jw 14 


A/21 


M 22 


MR 





A/31 


M 3 R 2° 








M R = 












+(-1 



a-b\ 















(1 








-l a A/i* 3 








A*2*2 


-l 6 A/ 2 *i 







"I'M*! 





A/11 A/{ 
A/ 22 MR 

A/Jy ±m. 









■ A/ 



3 2 







TMi* 







f Mn 








..rim 
Jw 14 




= 1 





M 22 


MR 





++- 




1 


M 3 R 2° 


±A/ 2 *2 









\M 4 ™ 








TM?i 



A/11 


A/12 


A/13 


.rim 
M 14 


A/21 


M 22 


i, rim 

Jw 23 


(1 


A/31 


MR 








M 4 » 












-(-l a - 6 ) 
























-l°Mi* 3 








M 2*2 


-l 6 A/ 2 *i 





-l"A/ 3 *i 


-l b M 2 *! 


A/*i 
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